Worm-like Polymer Loops and Fourier Knots 



O 
O 



o 

B 



X3 

o 

> 

in 
o 

00 

o 
o 



I 

o 
o 



X 
S3 



S.M. Rapaport'l', Y. Rabin''' and A. Yu. Grosberg-'- 
^Department of Physics, Bar-Ran University, Ramat-Gan 52900, Israel and 
^Department of Physics, University of Minnesota, 
116 Church Street SE, Minneapolis, MN 55455, USA 

Every smooth closed curve can be represented by a suitable Fourier sum. We show that the 
ensemble of curves generated by randomly chosen Fourier coefficients with amplitudes inversely 
proportional to spatial frequency (with a smooth exponential cutoff), can be accurately mapped on 
the physical ensemble of worm-like polymer loops. We find that measures of correlation on the scale 
of the entire loop yield a larger persistence length than that calculated from the tangent-tangent 
correlation function at small length scales. The conjecture that physical loops exhibit additional 
rigidity on scales comparable to the entire loop due to the contribution of twist rigidity, can be tested 
experimentally by determining the persistence length from the local curvature and comparing it with 
that obtained by measuring the radius of gyration of dsDNA plasmids. The topological properties of 
the ensemble randomly generated worm-like loops are shown to be similar to that of other polymer 
models. 



Classical polymer theories, from Debye to Flory to De 
Gennes |l[ and beyond, make extensive use of the analogy 
between conformation of a chain molecule and a Brow- 
nian random walk (BRW). In mathematically idealized 
form, such a random walk is thought of as a Wiener 
trajectory r(r) generated by the measure P {^{t)} oc 



exp 



const / {df/dr) dr 



where r is a parameter run- 



ning along the trajectory. The BRW model brings power- 
ful mathematical techniques, such as the diffusion equa- 
tion, to bear on polymers but, however fruitful, this 
method has its limitations. For example, it is funda- 
mentally unable to reproduce the fact of finite extensi- 
bility of polymer chains stretched by a strong forceQ- 
Another field in which the Wiener trajectory model fails 
miserably is the study of polymers with knots. Indeed, 
simulations of discrete polymer models, or random poly- 
gons with N steps, show that the probability of 
trivial knot configuration decays exponentially with N as 
Po{N) ~ exp {—N/N*), where N* defines the crossover 
from an unknotted to a knotted regime. It was argued 
that Wiener trajectory polymer models can not be used 
to calculate this probability 0| , since they correspond to 
the limit N — > oo, Z — > (Hs a step length); Nl"^ = const, 
in which the probability of a trivial knot vanishes. Notice 
that the contour length of Wiener trajectory diverges, 
Nl —^oo, so it is not surprising that these models do not 
exhibit finite extensibility. A better continuum model 
of polymers which can handle the constraints imposed 
both by finite extensibility and by the presence of knots, 
and can describe the elastic response of such objects, is 
that of a worm-like chain. However, while the statistical 
physics of linear worm-like polymers is well understood, 
little is known about the conformations and the topology 
of worm-like loops (WLL). Our goal in this article is to 
work out a method allowing computational generation of 
an ensemble of smoothly curved conformations of such 
polymers and to carry out the statistical analysis of their 



geometric and topological properties. 

The conformations of worm-like polymers 
are generated by the measure P{r(s)} oc 



exp 



-const 



(9^f/9s^)^ ds 



subject to the con- 



straint of non-extensibility, \dr/ds\ — 1. That simply 
means that conformations are Boltzmann weighted 
by the bending energy proportional to the squared 
curvature. In more general elastic models 0], the three 
generalized curvatures describing the object could be 
treated as independent Gaussian variables. However, 
this is true only for linear (with open ends) filaments 
and attempts to generalize these methods to the case of 
closed loop failed because the loop closure conditions 
introduce a non-local coupling between the curvatures 
that makes the problem practically intractable (except 
for the case of small fluctuations of a planar ring in 
which this coupling can be explicitly taken into account 
0). One natural way to generate the conformations of 
a closed loop is to expand each component (r) in a 
Fourier series: 
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27mr 
T 



Bt sin 



/ 2TmT\ 



(1) 



Here, r is a parameter along the curve r(r) (0 < r < T) 
and the summation goes up to some cutoff frequency 
f^max (if the series converges sufficiently rapidly, the cut- 
off can be replaced by infinity). Any conformation of 
a sufficiently smooth closed loop can be fully described 
by the set of Fourier coefficients {A^, B^} , giving rise to 
the concept of Fourier knots introduced in Notice 
that a Fourier knot is a much more general concept than 
the more traditional Lissajou knot, for which only one 
frequency is present for each coordinate direction. 

Using the above prescription one can generate an en- 
semble of loops of different shapes by taking the Fourier 
coefficients from some random distribution. However, 
even though all the loops have the same period T, their 
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contour lengths Lt = /q dT'\dr/dT'\ are different for 
each reahzation of the Fourier coefficients. Even more 
importantly, since, in general, \dr/dT\ ^ 1, loops gener- 
ated by 1^ do not obey the inextensibility condition and 
can not be used to model worm-like polymers. In order 
to generate an ensemble of different conformations of an 
inextensible loop of some well-defined contour length, one 
has to transform to a representation in which trajectories 
are parametrized by the arclength s 

s{t)= I dT'\dr/dT'\, with \dr/ds\ = l (2) 

and then bring all the generated conformations to the 
same contour length by a suitable afhne transformation 
of all lengths and coordinates (this transformation does 
not affect the topology since the latter is independent of 
the parametrization). Thus, we can replace Lt = s{T) 
by any standard length L, provided that we rescale all 
lengths using the afhne transformation: r yL/Lt- 
Therefore, in three steps (generation of random coeffi- 
cients and Fourier summation J^l, reparametrization |(2J), 
and afhne transformation s sL/Lt), of which the lat- 
ter two couple all Fourier harmonics together in a com- 
plex way, we obtain a statistical ensemble of smoothly 
bent conformations of an inextensible loop of length L. 
Let us now check if this ensemble is representative of 
physical conformations of a worm-like loop which pos- 
sesses some characteristic bending and (possibly) twist 
rigidities. 

The first step is to recognize that on length scales much 
larger than some microscopic cutoff, the conformations 
of a worm-like polymer are well represented by those of 
a BRW. For the latter, the amplitudes of the Fourier 
components can be readily shown to be of the form 
= X/n for n < n-a^a-K and zero otherwise, 
where A is a random number (say, uniformly distributed 
between —1 and 1; the choice of another symmetric in- 
terval is equivalent to rescaling the contour length, see 
below). As will be demonstrated in the following, short 
scale behavior characteristic of WLL can be obtained by 
replacing the abrupt cutoff rtmax by Fourier coefficients 
that decrease smoothly with a characteristic decay fre- 
quency no, 

{Aj„S:J^^^ = Ae^"/"V^, (3) 

(notice that the WLL and the BRW expressions for the 
coefficients coincide for n ^ ng ) . 

The next step is to notice that the conformations of 
linear worm-like polymer are governed by the persistence 
length, which marks the cross-over between straight line 
conformations at small length scales and random walk 
type behavior at large length scales m. We expect that 
on small scales, the conformations of linear worm-like 
polymers and WLL are quite similar and that for the 
latter this cross-over is represented in Fourier space in 



terms of riQ. Of course, the analogy is meaningful only 
as long as no ^ 1, when there is a sufficiently broad 
range of length scales between the persistence length and 
the contour length of the knot. The characteristic prop- 
erty of the worm-like chain model is exponential decay 
(as exp(— s/^i)) of the tangent-tangent correlation func- 
tion, which defines the persistence length li \ similarly, for 
WLL, on length scales sufficiently small compared to the 
entire loop, one expects that {{t{s) ■ t{0))) ~ exp(— s/^i), 
where (()) denotes averaging both over the contour of a 
given loop and over the ensemble of loops. This expec- 
tation is confirmed in Figure^, where the logarithm of 
the correlation function is plotted against the dimension- 
less arclength s/h. The choice li = 0.43/no allows us to 
superimpose data for different values of no in the range 
30 < no < 270 (the small shift of the exponent is the 
result of the finite discretization of the contour length s). 

At first sight, analogy with worm-like chain models 
of linear polymers suggests that on length scales much 
larger than li the conformations of the loops are those 
of BRW with step size (Kuhn length) given by twice the 
persistence length, 2li. However, this analogy is open 
to question because the large scale behavior of a loop is 
strongly affected by the loop closure constraint and it is 
not clear, a'priori, whether the same persistence length 
controls both the small and the large scale behavior. We 
therefore decided to apply our method to generate a rep- 
resentative ensemble of configurations of WLL and use it 
to compute the tangent-tangent correlation function and 
the mean squared distance between two wcll-scparatcd 
points s and s' (|s — s'\ 3> ^i) along the contour. 

A simple estimate shows that on length scales much 
larger than the persistence length, the tangent-tangent 
correlation function of a WLL approaches a constant 
negative value {{t{s) ■ t{0))) ~ —2I2/L where I2 is, in 
general, different from li (the negative value of the cor- 
relation follows from the fact that the tangent has to 
turn back on itself in order to come back to its initial 
direction upon traversing the contour of the loop). In 
Fig. this correlation function is plotted in the interval 
0.1 < s < 0.5 (here and in the following we take L = 1) 
for different values of no. Upon averaging the correlation 
functions over the oscillations, all the results for different 
values of no can be collapsed to a single horizontal line by 
dividing them by I2 — 0.575/no. We therefore establish 
that unlike the case of linear worm-like chains in which 
only a single persistence length exists, WLL have at least 
two distinct persistence lengths, with I1/I2 — 3/4. 

We now proceed to compute the mean squared 
distance between two points on the loop sep- 
arated by a contour distance s, ^(i?^(s))) = 

(1/i) {^/p^ |r (s' -f s) — r (s')l^ ds'^ Since we expect 

WLL to behave like BRW on scales much larger than 
some persistence length Z3, the probability distribution of 
-R^(s) can be easily written down assuming that both s ^ 
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FIG. 1: (Color online) (a) Universal semi-log plot of the short length scale behavior of the tangent-tangent correlation function 
vs the scaled distance, s/li\ the best fit exp[—[s/l\ —0.07)) is given by the solid line, (b) Large scale behavior of the correlation 
function, for different values of no (the height of the curves increases with no), (c) Scaled (by ^3) mean-squared distance 
between points on the loop. The BRW expression is plotted by the solid line. Different values of no are represented by black 
(no = 30), red (no = 50), green (no = 70), blue (no = 100), cyan (no = 150), magenta (no = 200), yellow (no = 220), dark 
yellow (no = 250) and navy (no = 270). 
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FIG. 2: (Color online). Log-log plot of persistence lengths 
hyhjs (black squares, green triangles and red circles re- 
spectively) vs no. The expressions 0.435/no and 0.57/no are 
shown as solid lines. The standard deviations are represented 
by error bars. 



I3 and L—s ^ ^3. Apart from a normalization factor, this 
probability is equal to the product of two Gaussian func- 
tions: exp[-i?V(2s4)]exp[-i?V(2(L-s)4)]. Aver- 
aging with this distribution yields the well known rela- 
tion (see, e.g., [HI), (i?^ (s))^^^ = 2I3S (1 - s/L) . This 
result indicates that the plot of {{R'^{s))) / I3L against 
a = s/L is universal, i.e., it is the unit height parabola 
4(t(1 — cr), independent of either persistence length I3 
or total contour length L. Figure QJ; shows that the 
data averaged over 1000 different configurations collapse 
quite accurately on the expected parabolic master curve. 
Furthermore, by looking at our computational data for 
((i?^ (s — i/2))) averaged over all pairs of diametrically 
opposite points on the loop (s = L/2 = 0.5) we were able 
to relate the persistence length to the cutoff no. Within 
the accuracy of our simulation ^3 coincides with I2 (see 
Fig. 



The conclusion that local and global statistical prop- 
erties of worm-like loops are characterized by two differ- 
ent persistence lengths, h and l2i= h), respectively, is 
a new result of this work. Although analytical theory 
of worm-like loops does not exist at present, we can offer 
some tentative considerations about the origin of the two 
persistence lengths, based on the study of small fluctua- 
tions of elastic rings |^ . For rings of zero linking number 
that possess both bending and twist persistence lengths, 
it was shown that while on small scales writhe fluctu- 
ations depend only on the bending persistence length 
(on such scales the ring is essentially straight and bend- 
ing and twist decouple), both bending and twist persis- 
tence lengths contribute to the writhe on larger scales, 
for which the global geometry of the closed ring becomes 
important. Extrapolating this result to the case of a 
strongly fluctuating loop considered in this work, sug- 
gests that while h represents the pure bending contribu- 
tion to the persistence length, I2 contains both the bend- 
ing and the twist contributions and is therefore larger 
than li 13]. 

We now turn to the topological properties of WLL and 
examine the dependence of the trivial knot probability Pq 
on the cutoff frequency no. Since no determines the Kuhn 
length defined as either 2/i or the 2I2, we can bring our 
knotting probability data to a form comparable to that 
for discrete polymer models, where trivial knot probabil- 
ity depends on the number of segments N. In our case, 
we generate Fourier knots as completely smooth curves, 
but we can define the number of effective segments as 
the ratio of the contour length to the Kuhn length. For 
long random- walk- like polymers, this definition coincides 
with the standard one accepted in polymer physics for 
the Kuhn segments lll| but there remains an ambiguity 
associated with the choice of A^i = L/2li or N2 = L/212- 
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FIG. 3: (Color online). Semi- log plot of the trivial knot prob- 
ability as a function of no. 

In order to resolve this ambiguity we will measure the 
probabihty to obtain a trivial knot as a function of the 
cutoff frequency uq. 

To address the topology of the loops computationally, 
we employ the knot analysis routine due to R. Lua 
which identifies knots by computing the Alexander poly- 
nomial invariant A(t) at one value of argument, t = —I, 
and Vassiliev invariants of degrees two, V2, and three, 
?;3. This set of invariants is widely considered as pow- 
erful enough for reliable identification of the trivial knot 
for all lengths achievable in practical computations. The 
details of this topological routine are described in . 

The data on the trivial knot probability are shown in 
Figure 13 As this figure indicates, the trivial knot prob- 
ability fits well to the exponential 

Po (no) ^ exp [-no/242] . (4) 

Notice that since this relation involves only the cutoff on 
the Fourier series, formula © can be re-interpreted in 
purely mathematical form, not involving any references 
to polymers, or, for that matter, to any physics. While 
there is no fundamental understanding of the origin of 
this large cutoff (rig = 242) at present, our formulation 
hints at the existence of a hitherto unexplored connection 
between Fourier analysis and topology of space curves 
and will hopefully stimulate new work on this fundamen- 
tal problem. 

In summary, we have shown that the choice of random 
Fourier coefficients with amplitudes that decay with fre- 
quency as exp(— n/no), generates a statistical ensem- 
ble of Fourier knots whose local properties coincide with 
those of worm-like polymers with persistence length that 
scales as I/uq. We found that even though WLL be- 
have on large scales as BRW as expected, the effective 
step size of this BRW is larger than that calculated from 
the local persistence length of the loop. This new pre- 
diction can be directly tested by experiments on double 
stranded DNA loops that can monitor both local (persis- 
tence length) and global (e.g., radius of gyration) prop- 



erties of the polymers (see, e.g., 0, 0, 0). We also 
demonstrated that similar to discrete models of polymers 
or random polygons, our Fourier knots exhibit exponen- 
tial decay of the unknotting probability with the number 
of effectively straight segments, or, equivalently, with the 
maximal spatial frequency included in Fourier expansion. 
The characteristic cross-over determined by this expo- 
nential decay represents a large number, which, although 
in the same ballpark as for other known models, remains 
an unexplained puzzle. 
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